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We calculate shear viscosity of the quark matter at finite temperature and density. If we assume 
that the quark interacts with the soft mode, which is a collective mode of quark-antiquark pair, the 
self energy of the quark is calculated by using the quasi-particle random phase approximation. It is 
shown that its imaginary part is large and the mean free path of the quark is short. With the use 
of the Kubo formula, the shear viscosity of quark matter becomes small. The Reynolds number of 
the quark matter is estimated to be about 3 ~ 30. As the temperature increases, the shear viscosity 
increases gradually for T > 200MeV. Moreover it is shown that the shear viscosity is not sensitive 
to the chemical potential. 
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I. INTRODUCTION 

The existence of quark matter or quark gluon plasma 
(QGP) has been predicted by Quantum Chromodynam- 
ics (QCD) which is the fundamental theory of hadronic 
world. It is supposed that the quark matter may ex- 
ist in the early stage of the universe or in the core of 
the compact stars. Recently such a new state of matter 
has been explored in the high-energy heavy ion collisions 
(RHIC) [I]- [!]• Now it is very important to study the var- 
ious physical properties of the quark matter in relation 
to these researches. 

Several years ago one of the present authors (M.I.) 
studied the specific heat and compressibility of the quark 
matter by using the Nambu-Jona-Lasinio (NJL) model 
[5]- These quantities are enhanced near the critical tem- 
perature. On the other hand the transport coefficients, 
shear viscosity and heat conductivity, are also important 
in relation to the high-energy heavy ion collisions. Re- 
cently it was discovered that the quark matter seems to 
be almost perfect fluid above the critical temperature. 
Therefore the calculation of the shear viscosity has been 
paid attention by many authors from the various view- 
points [6]-[26], [21]- [28] (lattice QCD), [29]- [31] (Black hole 
physics) . 

In this paper we try to calculate the shear viscosity 
of the quark matter by using the Kubo formula and the 
NJL model. We start to discuss the shear viscosity in 
terms of the Kubo formula [52"]-|35], which is expressed 
by the correlation function of the energy-momentum ten- 
sor. Suppose that the quark system in the QGP is a free 
gas. Then the shear viscosity becomes infinite because 
of the strong correlation between their current densities. 
On the other hand the RHIC data suggests very small 
viscosity. In order to clarify the physical origin of this 
small viscosity, we consider the shear viscosity of quark 
sector with the use of the NJL model, which does not 



include gluons explicitly. In the previous paper (referred 
to as I), the quasi-particle random phase approximation 
(RPA) was applied to this formula and it was shown that 
the shear viscosity is represented by the quadratic form 
of the quark spectral function |36| . Consequently it is 
the main purpose of this paper to calculate the spectral 
function within the NJL model and obtain the values of 
the shear viscosity under the various temperature and 
density. 

The spectral function is related to the self energy of 
the quark. In particular it depends on its imaginary part 
(or mean free path) strongly. What interaction does en- 
hance the self energy? Here we are reminded of the soft 
mode, which was proposed as a precursor of the chiral 
phase transition by Hatsuda and Kunihiro long time ago 
[57] , These modes grow strongly just above the critical 
temperature which is just the region we are interested 
in. We think that the collective mode couples strongly 
to the quark. We calculate the self energy of the quark 
due to the coupling to the soft modes. As is pointed out 
in [35], the spectral function of the soft mode decreases 
as the temperature increases. This means that the shear 
viscosity seems to be enhanced as the temperature in- 
creases. We will investigate the temperature dependence 
of the shear viscosity as well as the density dependence. 
Finally we calculate the Reynolds number of the quark 
matter in order to discuss the magnitude of the calcu- 
lated shear viscosity. 

In the next section, it is reviewed that the shear viscos- 
ity is expressed in terms of the quark spectral function, 
which is the main result of I. Then the spectral function 
is represented by the self energy of the quark (§.3). In 
§.4, the main part of this paper, we calculate the self en- 
ergy by using the coupling to the soft modes in the quark 
matter. The numerical calculations are carried out and 
their results are discussed in the final section (§.5). 



II. SHEAR VISCOSITY IN THE NJL MODEL 
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The Kubo formulas for the hydrodynamic transport 
coefficients are expressed in terms of the correlation func- 
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FIG. 1: The ring diagrams for the correlation function II with 
r = 1 or 175-7". 



tion of the energy- momentum tensor J xy (r, t). According 
to I, the formula for shear viscosity is given by 



rj = rj{uj = 0) 



duj 



Imn R (w) 



(1) 



where n R (w) is a retarded Green's function obtained by 
analytic continuation of the following (imaginary time) 
Green's function: n R (a->) = n(iw n )L 



n(iw„) 







dr(T T (J xy (r,T)J xy (0,0))) 



(2) 

Here cu n — 2imT (n — 0, ±1, ±2, • ■ • ) is the Matsubara 
frequency for bosons and T T means the (imaginary) time 
ordering operator. 

We take the NJL model for the quark matter in this 
paper |38j . Then the Lagrangian density is given by 

C = $(i7 • d - m)i/} + ,g[(^C) 2 + ($il5Tip) 2 ] , (3) 

where tp is the field operator for quarks and g is their 
coupling constant. We assume that the current quark 
mass is zero (m = 0) for convenience. The canonical 
energy-momentum tensor is read as 



(4) 



When this equation is substituted into Eq.(j2|, we have 
an infinite series for the correlation function II . 

According to I, we approximate the higher orders by 
the ring diagrams shown in Fig(I] (random phase approx- 
imation). As is shown in I, the higher order terms which 
consists of two loops or more vanish. Because the leftist 
(rightist) loop of each diagram contains odd number of 
7 p (/i = 0,1,2,3) matrices and its trace vanishes. Thus 
only the first term of Fig.l survives. The correlation 
function is written as 



IL(w n ) 



1 V 

ft 



(3 



d 3 p 



p x 2 Ti[j 2 G(p,iei+iuj n ) 



7 2 G(p,i£i) 



(5) 



where G(p, ie{) denotes the dressed quark propagator. If 
the free propagator is substituted instead of the dressed 
one, we have an infinite viscosity. So finite viscosity 
comes from the finite life time of the quark or finite mean 
free path in the present model. The spectral representa- 
tion for the dressed propagator is written as 



G Q/3 (p, isi) 



de p a p(p,e) 



2tt 



l£l 



(6) 



where e x = 27r(n + 1)T (n = 0, ±1, ±2, • • • ) is the Mat- 
subara frequency for fermions. Substituting this expres- 
sion into the correlation function Eq.([5]), we obtain the 
expression of rj in terms of quark spectral function as 
follows: 



dc 

2^ 



d 3 P 2 

a 0c 



dnp{s) 
de 



Tr[ 7 2 p(p,e)7 2 p(P^)] 
(7) 

Here Uy{s) = (1 + <C e ) 1 denotes the Fermi distribution 
function (/? = 1/T). The detailed derivation of this equa- 
tion was given in I. The factor dnp/de in the right-hand 
side implies that the shear viscosity is crucially affected 
by the distribution function near the Fermi surface. 

Next we calculate the spectral function of quarks. The 
spectral function is defined by 



p(p,e)^[G R (p, £ )-G A (p, £ ) 



(8) 



where the retarded and advanced Green's function are 
written as 

-rR, 



G*(p) 
G A ( P ) 



G(p,i£i : 

G(p,iei 



£ - 

■ £ ■ 



15) , 

16) . 



(9a) 
(9b) 



The full Green's function G is expressed in terms of the 
self energy S as follows 

G(j>)- l = GW{p)- l -Y,(p), (10) 

where G^°\p) denotes the free Green's function. Here 
we introduce the projection operators for positive and 
negative energy states, A± = ^(lT7°7'p)- According to 
H.AWeldon [39], the free Green's function G^ is written 
as a sum of two (positive and negative) parts. 

G(°)(p)- 1 = (p + +p)7 A + + (/ + -p)7°A_ , (11) 

where p° + = iei + ji (/i =chcmical potential). Similarly, 
the self energy can be written as 

S R (p) = E + (p)( 7 A + ) + E_(p)( 7 °A_), (12a) 
S A (p) = E;(p)( 7 A + ) + El(p)( 7 A_). (12b) 



By substituting Eq.(|ll| and (|12j into Eq.([10|, we obtain 
the expressions, 



G R (p) = 



G A ( P ) 



7 °A^ 
7 °A_ 



iS 



iS 



P° + 



7°A+ 
7°A + 



iS ' 
(13a) 

~5 ' 
(13b) 
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FIG. 2: The soft modes. 



Here we used the following formula, 



1 



a7°A + + &7°A_ a 



(14) 



where a and b are arbitrary variables. By using these ex- 
pressions, the spectral function p is rewritten as follows: 

p(p) = P + (P)7°A + + p_(p)7°A- , (15) 
where the coefficients p±(p) are defined by 



P±(p) = 



-2 ImSz, 



(^Tp-ReS T ) 2 + (ImE T ) 2 



(16) 



Thus the spectral function is written in terms of self 
energy. In other words, our calculation of the shear vis- 
cosity is reduced to that of the quark self energy. It is 
instructive to notice that if there is no interaction, the 
spectral function is simply given by 

p (0) (p,£) = 2ttS( P ° + -ph°A+ + 2tt5(p + +p) 1 °A^. (17) 

In this case, the viscosity r\ obviously diverges, which in- 
dicates that the quark interaction is essential for realizing 
finite viscosity. 



III. SOFT MODES IN THE QUARK MATTER 

In order to obtain finite shear viscosity, we have to 
know the self energy of quark £. As the largest contribu- 
tion to the self energy, we assume that the self energy is 
acquired by the coupling to the soft modes, which are pre- 
dicted in the scheme of the N JL model [37] , [ID] . The soft 
mode is a collective mode of quark-antiquark pair above 
T c as shown in Fig(2] The existence of such mode is one 
of the precursor phenomena of the chiral transition and 
it behaves as a bosonic quasi-particle. As the soft modes 
are carriers of chirality, isospin and energy momentum, 
it is likely that the modes affect various transport coeffi- 
cients including viscosity. 



The propagator represented by Fig(2]is written as 
D a (q) = D^\\ + X a D^ +X a D^X a D^ + ■■■) 

(18) 



D (0) 



1 -DWX a (q) ' 

with D(°) = —2g. The factor 2 in this equation is a 
symmetry factor of the interaction ( , ipT a ip) 2 . The polar- 
ization operator X a is defined by 

d 3 p 



(2tt) 3 

xTr[T a GW(p + q J «e J + iw ri )r a GW(p,«e J )] . (19) 

Here T° = 1, T a — ij 5 T a (a = 1,2,3) and the trace 
should be summed up over all Lorentz, color and flavor 
indices. Using the relation, T a ^ fl T a = 7^, the above po- 
larization operator and the propagator of the soft mode 
are independent of a: X a (q) = X(q) and D a (q) = D(q). 
The spectral representation of is read as 



G<°>(k,ie m ) = 



de p(°>(k,ie) 



(20) 



2tt ie m - e 

Substituting this representation into Eq.(19), we obtain 

d 3 p f de f de' 



2tt 



1 



(iei + iuj n - e)(iei - e') 
xTr[^(p + q,£)p (0) (p,e')]- (21) 
Using the summation formula of Matsubara frequencies, 
1 



TV. 

/L^l In 



(iei + iu> n - e){iei - e') 



1 



1 



tanh tanh — 1 



E 

2 e - e> - iuj n V 2T 2T J 

we obtain the expression of the polarization operator, 

. f d 3 p fde fde' 1 

X ^=J(2^j2nj2, 



> (22) 



1 L i. £ 

x - tanh tanh 

2 V 2T 2T 



e — £' — lUJ r , 
£ 



X 



Tr[/°)( P + q,e)p(°)( P , £ ')] • (23) 



Now we consider the self energy of quark acquired by 
the coupling to the soft modes shown in Fig[3] The self 
energy is written as 

d 3 A; 



s(p*)=X>E/(^ 



a m 



;r Q G(°)(k,ie m )r Q C( P -k,z £; -i£ m ) , (24) 



4 



I = 
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where we have introduced four- vector k± = (k, ±fc). If 
the z-axis is taken in parallel with p, the integral over k x 
(k y ) vanishes (k± ■ 7 — > ±£7° — k z j 3 ) and yields 



d 3 fc f dwImD R (p-k,w) 



coth h tanh 

2T 2T 



(2tt) 3 J 2tt 2 

fc - /A 7 - (p • k)p • 7 



w - 



fi - iei 



coth ^ _ tanh ^±iA 7° + (p-k)p- 7 
2T 2T J u-k-fi-iei 



(30) 



FIG. 3: Soft mode coupling. 

Here let us substitute the spectral representations of 
and D. The spectral representation of G^ is given in 
Eq.(20) and that of D is denoted by 

D(q,iuj n ) = - duj ; . (25) 



After these substitutions, the self energy Eq.(24l can be 
rewritten as, 

ImD R (p - fc) 



(fe m - e)(u; + «£,„ - «£;) 



(26) 



Here again we made use of the relation r Q 7 p r Q = 7^. 
Using another summation formula of Matsubara frequen- 
cies, 



1 



- e)(«e m + w - iej) 

1 1 / , w 
= coth — 

2 iej - uj - e V 2T 

we get the expression, 
d 3 fc 



tanh — j 
2TJ 



(27) 



S(p,iej) 



de 

(2tt) 3 7 27 
(°)(k,£)InxD R (p-k,a/) 



dw / , £ 
— tanh — — 
2vr V 2T 



This equation is a linear combination of the following 
form: 



E(p,fei) = S (p)7° + S 1 (p)p.7. 



(31) 



The first and second terms represent the scalar and vector 
parts of the self energy respectively. 



IV. APPROXIMATION METHOD 

Before calculating the shear viscosity, it is appropriate 
to carry out the following approximation without loss of 
the physical essence. In our calculations we have three- 
dimensional integration containing |p±q|. Then we take 
the following approximation: 

|p ± q| ~ p9(p -q)+ q6{q - p) = n(p, q) , 

where 8 denotes usual step function: 6(x) = 1 (x > 0) 
and 6{x) = Q(x < 0). This approximation is to replace 
the scalar product p • q by its average. From the above 
approximation, p • q vanishes because p • q = [(p + q) 2 — 
(p — q) 2 ]/4 ~ 0. With the use of this approximation, we 
can proceed our calculations furthermore. 



coth^) 



Soft mode 



where the factor 4 comes from the number of the soft 
modes. Here it is noted that the principal value for 
coth(o;/2T) should be taken in the integral. If we substi- 
tute the spectral function Eq.(17) given in the previous 
section into Eq.(26), the self energy is read as 

d 3 fc 
(2tt) 3 



(28) 



The polarization operator X(q) in Eq.(23l is repre- 
sented by multiple integral. Using our approximation, 
the trace of Eq.(23l can be written as follows: 



Tr[ i0 ( )(p + q ) £)p( )(p,£')] 



4tt 2 



E(p,ie,)=-4 



dto ImD R (p - k,w) 
2?r ~ 2k 



pn(p, q) 

x [S(e + n - n(p, q)) -5(e + ^+ n(p, q))] 
x [6(e' + n-p)- 5(e' + fj, + p)]. 



[(£ + M )(£' + m)-P 2 



(32) 



coth ■ 



2T 



coth tanh 

2T 



tanh 
k 



k — n 
2T 



k + ■ 7 



2T 



lj + k — /i — iei 
k- ■ 7 



fl - lEi 



Carrying out the angular integration of Eq.(23l, we get 

f c N f 



(29) X(q,iu m ) = 



2ir 2 



S(p, q)(n p + m p -l)p dp, (33) 
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where S(p, q) is defined by 
q-p i 1 



Sip, q) 



1 



1 \P-Q- lUr, 
'/ + !> I 1 



q \p+q- lljJm 
2 2 



2p — iuj m 2p + iuj r , 



p- q + iuj n 
1 

p + q + iuj„ 

eip-q). 



(q-p) 
(q-p) 

(34) 



In the right hand side, the n p = [e@( p ^ + 1] 1 and 
m p = [e^^ p+ ^ + denote the quark and antiquark 
distribution functions, respectively. Thus we can cal- 
culate the real part (ReX R ) and the imaginary part 
(ImJ^) of X R from this equation by using the relation, 
1/(2 + i5) — P/x — inSix). From Eq.(18l, the imaginary 
part of D R is read as 



ImD K (q,oj) 



-Ag 2 ImX R iq) 



(1 + 2gReX R iq)) 2 + 4.g 2 (ImA«(g)) 2 



This equation plays a fundamental role in the later dis- 
cussion as far as the soft mode is concerned. 



B. Self energy of quark 

Next we go on to the calculation of the self energy £(p) . 
From the above equation, we substitute p • k = into 
Eq.(30), the vector part of the self energy Si vanishes. 
Hence we get 



So(p)7° 



(36) 



This means that the self energy is independent of the di- 
rection of the momentum. The final expression for S (p) 
is 

S (p, l£; ) = -2 J ^3 J ^ImD K iK(p,k),e) 
k — /i x 



coth — + tanh ■ 



uj k + a 

coth tanh 

2T 



1 

to + k — /i — iei 
1 



2T J uj — k — ji — %E\ 



(37) 



The retarded function is obtained through analytic 
continuation (j£j — > s+id) of this equation. Its imaginary 
part is given as follows: 



Im£«(p,e) = — dkk 2 



coth 



uj — k + (i 
2T 



tanh J* ) lmD R in, oj - k + fi) 



, uj + k + 11 , k + u \ 
coth — - tanh ) ImZ? K (K, uj + k + fi) 



2T 



2T 



As for the calculation of the real part Re£^, which is 
the potential energy of quark, we make use of the disper- 
sion relation (Hilbcrt transformation) as follows: 



ReS^(p, £ ) 



ImX$ip,e') 



de', 



(39) 



(38) 



where we should take the principal value of the integral 
in the right hand side. Now that we have obtained the 
self energy S, we can calculate the shear viscosity using 
Eq.(7). 



V. NUMERICAL RESULTS 

We are now ready to evaluate the shear viscosity at 
various temperature and density. The parameters of the 
NJL model are taken to be the same values as those in the 
reference [35]: the coupling constant, g — 5.2 MeV and 
the cut off parameter, A — 631 MeV. The quark mass 
m is assumed to be zero so that the chiral symmetry 
holds exactly. First we restrict ourselves to the case of 
the chemical potential /i — 0. Note that the RHIC data 
covers fi ~ 10 MeV. Under these parameters, the critical 
temperature T c is 161 MeV, which is derived from the 
vanishing condition of the inverse propagator of the soft 
mode: D" 1 (p = 0) = [35J- 



A. Soft mode 

To begin with, let us consider the soft mode in the 
quark matter. Its spectral functions ImZ? R (p = 0, uj) 
are calculated from Eq.(35) and the results are shown in 
Fig.4 at various temperatures: T = 180, 210, 260 MeV. 
For T very close to T c , the soft mode becomes a sharp res- 
onance and acquires the nature of a good quasi-particle. 
As the temperature increases, the resonance width is 
spread out and its strength is weakened at 170 MeV < 
T < 260 MeV. As is pointed out by Hatsuda and Ku- 
nihiro, this peak near the critical temperature reflects 
a precursor of the chiral phase transition. Therefore 
its existence implies that the system is in the restored 
chiral symmetry phase or QGP phase. In Fig. 5, we 
present the spectral functions with various momenta at 
T = 170 MeV. The position of the peak does not depend 
on the momentum p but the strength depends on p or 
the temperature T. Hence the soft mode is determined 
by the temperature. It will be shown later that the soft 
mode is not sensitive to the chemical potential. 



B. Self energy of quark 

Next let us calculate the self energy of the quark, which 
is described by the coupling to the soft modes as dis- 
cussed in the previous section. The real and imaginary 
parts of E do not depend on the momentum p sensitively 
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FIG. 4: The spectral function of the soft mode at T = 
170, 210, 260 MeV as a function of the frequency of the soft 
mode. 



FIG. 6: The spectral function of the quark at the momenta, 
p = 0, 200, 400 MeV (T = 180 MeV), which is generated by 
the coupling to the soft mode. 
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FIG. 5: The spectral function of the soft mode at the mo- 
menta of quark, p — 0, 110, 150 MeV as a function of the fre- 
quency of the soft mode (T = 170 MeV). 



as well as the soft mode. We show the spectral function 
of quark in Fig.6. They have a broad and small peak 
except for the case of p = 100 MeV. The position of the 
broad peak moves with the momentum such that the de- 
nominator in Eq.(16) is minimized. The sharp peak near 
e <~ at p — 100 MeV seems to come from the over- 
lap between two spectral functions, and lmD R . We 
see that the strength of the spectral functions are sup- 
pressed due to the integration of the momentum k. This 
behavior lowers the shear viscosity as shown in the next 
subsection. 



C. Shear viscosity 

Now we can calculate the shear viscosity by using 
Eq.(7). The results are shown in Fig. 7 as a function of the 
temperature above the critical one, T > T c = 161MeV. 
Similarly its ratio to the entropy is drawn in Fig. 8. The 
solid lines in these figures represent guidelines for the 
data. These values are almost constant but increase 
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FIG. 7: The shear viscosity as a function of the temperature 
T (n = 0). The dashed line denotes for result by the particle 
mode explained in the text. 
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FIG. 8: The ratio of the shear viscosity to the entropy as 
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denotes for result by the particle mode explained in the text. 
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FIG. 9: The spectral function of the quark at the momenta, 
p = 0,200, 400 MeV (T = 180 MeV), which is generated by 
the coupling to the particle mode. 



gradually at T > 200 MeV. The smallness of the shear 
viscosity is understood by the smallness of the quark 
spectral function as shown in Fig. 6. This comes from 
the strong coupling to the soft modes, which gives rise to 
the large self energy, that is the short mean free path. 
As the temperature increases, the soft mode becomes 
weak and its peak moves to the high frequency (Fig. 4). 
Then the quark spectral function is enhanced. This is 
the reason that the shear viscosity increases gradually at 
T = 280 MeV. 

For the comparison, we also calculated the shear vis- 
cosity by the coupling to the particle mode. This mode 
represents that the quark collides with another particle 
instead of the soft mode. The Feynman diagram of this 
process corresponds to the second diagram of the Fig. 3. 
Hence the spectral function of this mode is given by 



FIG. 10: The shear viscosity as a function of the chemical 
potential p (T = 180 MeV). 



ImD R (q) = 4(2g) 2 ImX R (q). 



(40) 



The calculated result is denoted by the dashed line in 
Fig. 7 and 8. We see that the viscosity by the particle 
mode is one-order larger than those by the soft mode. 
This means that the particle mode gives rise to the weak 
correlation between quarks and makes the quark mat- 
ter Fermi gas-like. This property is seen in the spectral 
function of the particle mode as shown in Fig. 9. 



D. Finite density 

Next let us discuss the density or chemical potential 
dependence of the shear viscosity. The shear viscosity 
for < fi < 100 MeV is shown in Fig. 10. The viscosity is 
increasing very slowly with the chemical potential. This 
dependence is understood by noting that the soft mode is 
not changed for the chemical potential as shown in Fig. 11 
at T= 170 MeV. 
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FIG. 11: The spectral function of the soft mode at the chem- 
ical potential, p = 0,50, 100 MeV as a function of the fre- 
quency of the soft mode (T = 170 MeV). 



E. Reynolds number 

Finally we will discuss how small the shear viscosity 
obtained in this paper is. To this end, let us introduce 
(dimcnsionlcss) Reynolds number which is defined by 

inertia term I ■ u 



R = 



viscosity term v/(pm) 



(41) 



In the numerator, I and u denote typical size and veloc- 
ity of the system, respectively. We set I = lOfm and 
u — 1 — light velocity for convenience. It is difficult 
to define "mass density p m " appeared in the denomina- 
tor. Here we have denned it by the qusi-particle mass m 
times number density of the quark and antiquark. If we 
set m — 10 <~ 100 MeV, the Reynolds number becomes 



R = 3 - 30. (for soft mode) 



(42) 



On the other hand, R = 0.3 ~ 3 in the case of the cou- 
pling to the particle mode. Taking into account that the 
shear viscosity is negligible under the large R, the evalu- 
ated R implies that our quark matter is close to a perfect 
liquid due to the strong coupling to the soft mode. 
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F. Summary 

In conclusion, we have evaluated the shear viscosity 
of the quark matter (QGP) by using the Kubo formula 
and the NJL model. The self energy of quark is calcu- 
lated by the coupling to the soft mode which is a collec- 
tive mode composed of the quark-antiquark pair. This 
strong correlation of quarks lowers the shear viscosity 
of the quark matter and the quark matter is nearly a 
perfect fluid. As the temperature increases, the shear 
viscosity increases gradually because the strength of the 
soft modes decreases with the temperature. On the other 
hand it is shown that the shear viscosity is not sensitively 



depend on the chemical potential. 
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